Measurement of the A and A particles 
in Au+Au Collisions at y^s^ = 130 GeV 
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We present results on the measurement of A and A production in Au+Au collisions at 
^s„„ — 130 GeV with the PHENIX detector at RHIC. The transverse momentum spectra were 
measured for minimum-bias and for the 5% most central events. The A/A ratios are constant as a 
function of pr and the number of participants. The measured net A density is significantly larger 
than predicted by models based on hadronic strings (e.g. HIJING) but in approximate agreement 
with models which include the gluon junction mechanism. 

PACS numbers; 25.75.Dw 



In this paper we report on the measurement by the 
PHENIX experiment at the Relativistic Heavy Ion Col- 
lider (RHIC) of the production of A and A particles and 
the ratio A/ A as a function of transverse momentum pT 
and centrality in Au+Au collisions at ^/snn — 130 GeV. 
The production of strange baryons and of strange par- 
ticles in general has been extensively studied in heavy 
ion collisions at the Alternating Gradient Synchrotron 
(AGS) and at the Super Proton Synchrotron (SPS) @] 
to investigate the flavor composition of nuclear matter at 
high density and temperature. Furthermore, antibaryon- 
to-baryon ratios, or alternatively, a net baryon num- 
ber such as (A — A) at midrapidity provide insight into 
the baryon transport mechanism in these collisions P-p| . 



The systematic study of baryon stopping (transport of 
baryon number in rapidity space) and hyperon pro- 
duction in proton-proton, proton-nucleus and nucleus- 
nucleus collisions has been done over the past decade at 
the AGS |,| and CERN SPS |-|l§. The results have 
shown a high degree of baryon stopping and enhanced 
hyperon production in heavy nucleus-nucleus collisions. 
Clearly the measurement of these quantities at RHIC, at 
the highest energies so far available in the laboratory, is 
of great importance for our understanding of the sources 
of these processes. 

The results reported here were obtained using the west 
arm of the PHENIX spectrometer which covers an 
angular range of = 7r/4 (during its first year of run- 
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ning) and a pseudorapidity range of \r]\ < 0.35. The 
detectors used were the drift chamber (DC), a set of 
multiwire proportional chambers with pixel-pad read- 
out (PCI) and a lead-scintillator electromagnetic 
calorimeter (EMCal) . Signals from two sets of beam- 
beam counters (BBC) and two zero-degree calorimeters 
(ZDC) provided a trigger sensitive to 92 ± 4% of the 6.8 
barn total Au+Au cross section [^. The centrality se- 
lection was done using the correlation between the analog 
responses of the ZDC and BBC @. 

The present analysis is based on 1.3 M minimum-bias 
events with a vertex position of \z\ < 20 cm. To recon- 
struct the A and A particles, their weak decays A pn^ 
and A pTT"*" are used. The tracks of the charged par- 
ticles from the decay of A and A are reconstructed using 
DC and PCI, and their momentum is determined by the 
DC with a resolution of ,5p/p ~ 0.6% ® 3.6% p (GeV/c). 
Although these tracks do not point back to the primary 
vertex position their momentum is calculated assuming 
that they come from the interaction point, hence there is 
in general a shift in the momentum of the reconstructed 
A and A particles. A Monte Carlo (MC) study shows that 
the difference is of the order 1-2% over the measured mo- 
mentum range, within the measured momentum resolu- 
tion, and thus neglected in the present study. The abso- 
lute momentum scale is known to better than 2% . In 
order to reduce background the tracks are confirmed by 
requiring a matching hit in the EMCal within ±3cr. For 
the particle identification (PID) of the daughter charged 
particles the time-of-flight signal of the EMCal with a 
timing resolution of ^ 700 ps is used. Using the momen- 
tum measured by the DC and the flight-time, the particle 
mass is calculated and a 2a momentum-dependent cut is 
applied to the mass-squared distribution to identify pro- 
tons, antiprotons and pions. An upper momentum cut 
of 0.6 GeV/c and 1.4 GeV/c for pions and protons re- 
spectively provides clean particle separation. Then each 
proton is combined with each pion in the same event 
and the invariant mass is calculated. If A or A are pro- 
duced a peak appears in the mass distribution on top of 
a background from random combinations of particles. In 
order to determine the number of A's and A's from such 
a distribution an estimation of the background is essen- 
tial. For this, the mass distribution with combinations 
of protons and pions from different events with the same 
centrality class (event mixing) is used. In order to reduce 
the combinatorial background the decay proton energy is 
required to be within £;™™ < Ep < E"™"^ where E^'"- 
and E"™"^ are calculated from the two-body decay kine- 
matics in the A center-of-mass system. A similar cut is 
used for the pions. This cut gives an improvement of the 
signal-to-background ratio by a factor of two and results 
in the final value of S/B = 1/2 for both A and A. We 
obtain ^ 12000 A and ~ 9000 A particles in the mass 
range 1.05 < m^^r < 1-20 GeV/c^ with mass resolution 
5m/ m ~ 2% obtained from a Gaussian fit. 



Fig. |I| shows the invariant mass spectra for the A 
p'K~ and A piT~^ . The results represent the primary A 
and A and contributions from the feed-down from heavier 
hyperons (mainly E" and S). The reconstructed number 
of A and A particles is corrected for the acceptance, pion 
decay-in-flight, momentum resolution and reconstruc- 
tion efficiency. For this, single-particle MC events were 
generated over the full azimuth (j> {0 < (j) < 27r) and one 
unit of rapidity y (-0.5 < y < 0.5). The simulated par- 
ticles were passed through the entire PHENIX GEANT 
p6| simulation. The correction function is defined as 
the ratio of the input (generated) transverse momentum 
distribution to the pt distribution of the particles recon- 
structed in the spectrometer. However this correction 
does not take into account the decreasing track recon- 
struction efficiency due to the high multiplicity environ- 
ment in more central events. A well established method 
to obtain this efficiency drop is the embedding procedure 
Ip^ . We use single-particle MC tracks embedded into 
real events and analyze the merged events with the same 
analysis code as used for the reconstruction of the data 
set. The track reconstruction efficiency decreases from 
90% for minimum-bias to 70% for central events. The 
S° and S hyperons decay to A with a branching ratio 
of essentially 100%. We have verified, using HIJING, 
that the kinematic properties {px and y distributions) of 
primary A and those produced by and S decay are 
the same (within a few percent). We conclude that the 
correction function determined for primary A is valid for 
all A. Since there are no reliable data available for the 
yields of those hyperons at RHIC energies, we cannot 
quantitatively state the contributions from those heavier 
hyperons to our A production. Therefore, our data in- 
clude the primordial A and A, as well as the feed-down 
from the heavier hyperons. 

Using the correction function determined from single- 
particle MC and multiplicity-dependent track recon- 
struction efficiency derived from the embedding proce- 
dure we correct the transverse momentum distributions 
of A and A. The invariant yields as a function of the 
transverse momentum px for A and A are shown in Fig. ^ 
for minimum-bias events (circles). Over the measured 
transverse momentum range (0.4 < pT < 1.8 GeV/c) 
both the A and A spectra can be described by a distri- 
bution of the form dN/dpT^ cx p^e^^^^"^ as shown in 
Fig. m by the solid line. The total yield dN/ dy is obtained 
by integrating the functional fit from zero to infinity. In 
the same Fig. ^ we show the invariant yield for the 5% 
most central events (squares). In Table | the total yields 
dN / dy together with the average transverse momentum 
<Pt> and the slope parameters T are listed for both 
minimum-bias and for the 5% most central events. 

There are several sources which contribute to the sys- 
tematic uncertainties of the total yield dN/ dy and aver- 
age transverse momentum <pt> determination. One of 
the sources, the uncertainty in the correction function de- 
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termination, is found to be 13%. A second contribution is 
due to the fitting function used for the extrapolation be- 
yond the measured pr range: the fraction of the extrap- 
olated yield is 29 ± 9 %. There is also an additional con- 
tribution to the systematic error which originates from 
the combinatorial background subtraction which is 3%. 
Combining the above uncertainties in quadrature gives a 
total systematic error in the yield of 16% (see Table |). 

The A/ A yield ratio determined in the mass win- 
dow defined above versus px is shown in Fig. || (top 
panel). The ratio is constant over the whole px range 
0.4 < pt < 1.8 GeV/c. There is also no significant 
variation of the A/ A ratio as a function of the number 
of participants (bottom panel) which is calculated us- 
ing a Glauber model together with a simulation of the 
ZDC and BBC responses Q. The average A/A ratio 
is 0.75 ± 0.09 (stat). Both the pT-dependence and the 
integral A/ A ratio are consistent with statistical thermal 
model calculations [0 at RHIC energies. 

The present measurement of the total yield of A and 
A enables us to take the previously reported inclusive p 
and p spectra |l^ , and construct a feed-down correction 
for A decays. As the S yield has not been measured, we 
do not inlcude feed-down from decays, but this is ex- 
pected to be < 5%, based on HIJING calculations. The 
feed-down corrections were done bin by bin on the proton 
(antiproton) px spectrum by the following method: 

(0.1) 

where dN^ / dydpx is the total yield of the primary pro- 
tons, dN™' / dydpx - the total yield of the measured pro- 
tons, dN^ / dydpx - the total yield of the measured A, BR 
- branching ratio of the A decay A — > p7r~ , i is the px bin 
number, N{,ins is the number of bins, w{j, i) is the frac- 
tion of protons from A decay from bin number j which 
fall into the proton bin number i. These fractions were 
extracted from MC. The feed-down corrected proton and 
antiproton px spectra are shown in Fig. ^. We calculated 
the total yield, dN/dy, for protons and antiprotons by fit- 
ting them to the same distribution as used for the A and 
integrating from zero to infinity. The results for the total 
yields dN/dy, the average transverse momentum <px> 
and the slope parameters T for minimum-bias and for 
the 5% most central events are also listed in Table |. The 
measured A/p and A/p ratios after feed-down corrections 
are found to be 0.89 ± 0.07 (stat) and 0.95 ± 0.09 (stat). 

The net baryon numbers are indicative of the baryon 
transport mechanism in relativistic heavy ion (RHI) col- 
lisions. In Table ^ we compare our results for minimum- 
bias and for the 5% most central events with the predic- 
tions of the HIJING |l^ model which assumes that the 
primary mechanism for baryon transport in RHI is due to 



quark-diquark hadronic strings. Vance et al. |^ imple- 
mented a non-perturbative gluon-junction mechanism in 
a new version of HIJING (called HIJING/B) to explain 
the enhanced baryon stopping at CERN SPS energies. 
The predictions of this model are also shown in Table ||. 
For a valid comparison with the experimental data, the 
HIJING and HIJING/B results for A and A also include 
the feed-down from the heavier hyperons, and the results 
for protons and antiprotons the feed-down from . The 
Table shows a clear difference (by a factor of ~ 4) between 
HIJING and HIJING/B for the net A number, with the 
HIJING/B predictions in much better agreement, in par- 
ticular for the minimum-bias data, where the experimen- 
tal errors are relatively small. The difference between the 
two models for the net proton number is less obvious and 
although the data seem to favor HIJING/B the present 
accuracy does not allow for a clear preference between 
the two models. 

In conclusion, we have measured the transverse mo- 
mentum spectra of the of A and A particles in the px 
range 0.4 < px < 1.8 GeV/c in minimum-bias and the 
5% most central Au4-Au collisions at .^s„„ = 130 GeV 
at RHIC with the PHENIX experiment. The absolute 
yields of dN/dy, at mid-rapidity, of A hyperons are de- 
termined by extrapolating to all values of px- The av- 
erage A/A ratio is found to be 0.75 ± 0.09. The ratio 
is constant over the whole px range and there is also no 
significant variation of the A/ A ratio as a function of the 
number of participants. Using the measured A and A 
yields, the p and p yields corrected for feed-down from 
A decays are determined. The A/p and A/p ratios after 
feed-down corrections are found to be 0.89 ± 0.07 and 
0.95 ± 0.09. The measured net A number is substan- 
tially larger than predicted by HIJING as already seen 
at CERN SPS and which may indicate enhanced baryon 
stopping at RHIC energies. The newly available high 
statistics data at ^s„„ = 200 GeV will allow us to fur- 
ther study baryon transport and strangeness production 
in RHI collisions. 
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FIG. 2. Transverse momentum spectra of A and A for min- 
imum-bias and for the 5% most central events. For clarity 
of presentation the data points for minimum-bias are scaled 
down by a factor of ten. 
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FIG. 1. Invariant mass spectra of pTi~ (solid squares) and 
pTT^ (open squares) pairs. The histograms show the back- 
ground for pTx~ (solid) and pix^ (dashed). For clarity the 
p-K^ data are scaled up by a factor of two. The insert shows 
the background-subtracted spectra of A (solid circles) and A 
(open circles) . The lines are Gaussian fits to the A (solid line) 
and A (dashed line) mass spectra. 
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TABLE I. Inclusive A and feed-down corrected proton 
yields, average transverse momentum (in GeV/c) and slope 
parameters (in MoV/c) for minimum-bias (MB) and for the 
5% most central events. The errors listed are statistical. The 
systematic errors are 17% for the protons and 16% for the A. 



p 




















A 


A 


P 


P 






dN/dy (MB) 


4.8 ± 0.3 


3.5 ± 0.3 


5.4 ± 0.2 


3.7 ± 0.2 






(5%) 


17.3 ± 1.8 


12.7 ± 1.8 


19.3 ± 0.6 


13.7 ± 0.7 



<PT> (MB) 1.06 ± 0.08 1.10 ± 0.12 0.88 ± 0.04 0.91 ± 0.06 
(5%) 1.15 ± 0.15 1.14 ± 0.21 0.95 ± 0.07 1.04 ± 0.12 



(MB) 
(5%) 



355 ± 11 
384 ± 16 



366 ± 13 
380 ± 19 



292 ± 21 
319 ± 31 



304 ± 23 
327 ± 34 



FIG. 4. Spectra of inclusive A (A) and feed-down corrected 
protons (antiprotons) for minimum-bias events. The data 
points for antiprotons and A are scaled down by a factor of 
two. 



TABLE II. Total measured and predicted baryon yields for 
minimum-bias (MB) and for the 5% most central events. The 

errors listed are statistical. The proton yields are corrected 
for feed-down. The systematic errors are 24% for the protons 
and 23% for the A. 



Net baryon number PHENIX HIJING HIJING/B 



(A -A) 


(MB) 


1.3 ± 0.4 


0.2 


0.8 


(p-p) 


(MB) 


1.7 ± 0.3 


1.1 


1.7 


(A -A) 


(5%) 


4.6 ± 2.5 


0.8 


3.2 


(p-p) 


(5%) 


5.6 ± 0.9 


4.7 


7.1 
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